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ABSTRACT 

We explore the implications of models of the Hyades moving group in which it has a 
resonant origin, for regions of the Galaxy beyond the Solar neighbourhood. We show 
that while models associated with different resonances can produce nearly identical 
substructure in the local velocity distribution, the velocity distribution away from the 
Solar neighbourhood has different properties for different models. In particular there 
is a variation between different models of where in Galactocentric radius the observed 
Hyades signal in velocity space is strongest, at a given Galactic azimuth. We note, 
however, that the uncertainties in currently available data, primarily due to uncertain 
distances to stars, hide these signatures rather effectively, meaning we are not yet able 
to determine which resonance is the cause of the Hyades. 
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1 INTRODUCTION 

It is well established that the velocity distribution in the 
Solar neighbourhood is far from smooth. The existence of 
"moving groups" of stars with similar velocities in the Solar 
neighbourhood has been recognised for over a hundred years 
(Kapteyn 1905; Eggen 1996, and references therein), but it 
was generally believed that these groups are composed of 
stars from dissolved clusters, which retain similar velocities 
' (at a given point in the Galaxy) while spreading in space. 
Using observations by the Hipparcos satellite (ESA 
1997), Dehnen (1998) showed that these moving groups 
dominate the distribution of blue (B — V < 0.4) stars, 
and contain a substantial proportion of the redder stars in 
the Solar Neighbourhood (beyond Parengo's discontinuity, 
B — V > 0.61), with the same moving groups being seen in 
the velocity distributions of the different subsamples of Solar 
Neighbourhood stars. This suggested that a single moving 
group contains stars of many different ages, a result con- 
firmed by Famaey et al. (2005), who added to the Hippar- 
cos data by determining radial velocities for a sample of 
Hipparcos stars, allowing them to directly associate a given 
star with a given moving group. Comparison of these stars 
to isochrones in the Hertzsprung-Russell diagram then re- 
vealed the wide range of ages in each moving group. This is 
a very strong argument that moving groups must have dy- 
namical origins. More recently, Pompeia et al. (2011) used 
chemical "tagging" to show that a large fraction of the stars 
in the Hyades moving group can not have originated from 
the Hyades cluster. 
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The dynamical processes that shape the velocity dis- 
tribution of the Solar neighbourhood in to moving groups 
is still unclear. Kalnajs (1991) suggested that the Hyades 
and Sirius moving groups could be associated with an outer 
Lindblad resonance with the Galactic bar, while Dehnen 
(1999; 2000) used test-particle simulations to show how this 
resonance with the bar could instead produce the Hercules 
moving group (which is at velocities further from the circu- 
lar velocity than the Hyades or Sirius groups, and was not 
fully recognised when Kalnajs was writing). Fux (2001) came 
to a similar conclusion, and more recently Minchev et al. 
(2010) argued that the moving groups at velocities near to 
the circular velocity may also be due to the influence of the 
bar's OLR. Other work, beginning with De Simone, Wu, & 
Tremaine (2004) , has looked at the possible influence of spi- 
ral structure, and shown that it is capable of producing a 
variety of substructures in velocity space, including those ob- 
served in the Solar Neighbourhood (e.g Quillen & Minchev 
2005; Antoja et al. 2011; Pompeia et al. 2011). Other studies 
have sought to explain at least some of the moving groups as 
the response of the disc to the accretion of a satellite galaxy 
(Quillen et al. 2009; Minchev et al. 2009) or as the debris of 
an accreted satellite (Helmi et al. 2006, though this is only 
plausible for the higher velocity moving groups) . All of these 
works used some form of test particle integration or N-body 

1 In the discussion section of these conference proceedings, 
Kalnajs (questioned by Lindblad) concedes that the same struc- 
ture in the velocity space at a single point could be produced by 
an inner Lindblad resonance with some unknown perturbation, 
but notes "I could tell the difference if I knew the stellar velocity 
distribution in a large enough region around the sun. But I will 
never know that." 
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simulation to model the effects of these perturbations and 
their resonances. 

More recent work has used angle-action coordinates to 
show that the Hyades moving group 2 is associated with trap- 
ping at a Lindblad resonance (Sellwood 2010). However, it 
is not clear which resonance is responsible, as very different 
resonances are associated with very similar structures in the 
velocity distribution of the Solar neighbourhood (McMillan 
2011b, henceforth Mil). It is this approach that we pursue 
in this study, because it allows us to create easily tuneable 
models which make no assumptions about the nature of the 
perturber. 

Recently it has become feasible to examine the veloc- 
ity distribution of stars some distance away from the So- 
lar neighbourhood. Notably (Antoja et al. 2012, henceforth 
A12) used observations by the RAdial Velocity Experiment 
(RAVE: Steinmetz et al. 2006) to investigate the veloc- 
ity distribution in fields up to ~ 1.5 kpc from the Sun, in 
the anti-rotation direction. In a spatial bin located around 
the Sun's Galactocentric radius, but at least 700 pc away 
in the Galactic azimuthal direction, they found a structure 
in velocity space at a similar position to that occupied by 
the Hyades in the Solar neighbourhood, as well as other 
overdensities at other points in velocity space. In other spa- 
tial bins which lie at similar Galactic azimuths, but at least 
300 pc inside or outside the Solar radius, they find no such 
structure at the velocity of the Hyades, but still find other 
substructure in velocity space. 

In this paper we extend the models used by Mil to 
explore the implications of the expected velocity distribu- 
tions away from the Solar neighbourhood. In Section 2 we 
discuss the nature of the resonances considered and the use 
of angle-action coordinates in describing them. In Section 3 
we give the numerical details of the Hyades models that we 
test, before describing their properties away from the Solar 
neighbourhood in Section 4. 



2 ANGLE-ACTION COORDINATES AND 
LINDBLAD RESONANCES 

Three actions J; and three conjugate angle coordinates 9i 
provide exceptionally convenient coordinates for objects or- 
biting in a stationary or slowly evolving gravitational poten- 
tial (and natural coordinates for perturbation theory). The 
actions are conserved quantities and the angles increase lin- 
early with time, 6i(t) = #;(0) + fli(J)t, where f2; is a fre- 
quency. This means that J can be thought of as labeling an 
orbit, and as describing a point on that orbit. The usual 
phase space coordinates x,v are 27r-periodic in each angle 
coordinate For a phase-mixed distribution function (df) 
the distribution of stars is uniform in angle and therefore 
the DF / = /(J). 

The conversion between the coordinate systems 6, J 
and x,v is only known analytically for a few types of 



2 Throughout this paper we will refer to dynamical structure 
associated with the Hyades moving group as "the Hyades" in the 
interests of brevity. This should not be interpreted as meaning 
the Hyades star cluster, which is part of the moving group, but 
certainly doesn't extend beyond the Solar neighbourhood! 



gravitational potential. However we can use the "torus fit- 
ting" method (e.g. McMillan & Binney 2008, and references 
therein) to approximately find this conversion in an axisym- 
metric Galactic potential. This is done numerically, on an 
orbit-by-orbit basis - i.e. for a given J the torus fitting al- 
gorithm finds an expression giving the values of x, v for any 
0, but a new torus fit has to be performed to know this 
relationship (x, v for any 0) for a different J. 

In this study we work under the approximation that we 
can ignore the influence of the non-axisymmetric structure 
of the disc except for their trapping stars around resonances. 
This means that we can use the angles and actions calculated 
in the underlying axisymmetric potential, but that the df 
is reshaped by resonant trapping, and is not simply f(J) 
but also depends on the angles. This reshaping produces 
the Hyades. 

Sellwood (2010) explored the distribution of the Solar 
neighbourhood stars in angle-action coordinates (under sim- 
ilar approximations), using data from the Geneva Copen- 
hagen Survey (GCS, Holmberg, Nordstrom, & Andersen 
2009). The distribution of stars in action showed a clear indi- 
cation that the Hyades moving group was associated with a 
resonance between the radial and azimuthal orbital frequen- 
cies of the stars (Q r and Q4,, both of which are functions of 
J) associated with a perturbation whose pattern speed fl p 
satisfies 

lQ r (J) + mn^(J) = mQ p . (1) 

I defines the type of resonance, with I = +1 for an outer 
Lindblad resonance (OLR), I = — 1 for an inner Lindblad 
resonance (ILR), I — at corotation, and / = ±| for ultra- 
harmonic resonances. (Throughout this paper we will refer 
to a Lindblad resonance for any m as an m:l ILR or OLR.) 3 

The condition given by equation 1 is satisfied along 
planes in action space, and to a good approximation these 
can be thought of as a relationship between the radial ac- 
tion ,J r and the azimuthal action (which is the angular 
momentum about the symmetry axis of an axisymmetric po- 
tential), with the third action J z having a negligible effect 
on the relevant frequencies. It therefore makes sense to think 
about resonant lines in the J r -.J^ plane and ignore changes 
with J z . For the values of that are relevant for disc stars 
near the Solar neighbourhood, changing £l p shifts the zero 
intercepts of these lines (the value of for J r — 0, i.e. the 
angular momentum of the relevant circular orbit), but has a 
negligible effect of the slope of the line in the — J r plane. 

The resonant lines for different ILRs or OLRs that go 
through the Solar neighbourhood have very similar slopes, 
and the stars associated with the Hyades clearly lie around 
one of these, but it is impossible to tell from the distribution 
in J which values of I and m (and thus il p ) define the reso- 
nance (see Figure 1). It is also worth noting that the slopes 
of the lines in J associated with the various resonances can 
be seriously affected by relatively minor changes in the as- 
sumed potential. The two panels of Figure 1 illustrate this 
by showing the difference (most notably for the 2:1 ILR res- 
onance) between the lines assuming a logarithmic potential 



3 Note that in these models there is nothing to distinguish the 
effects of a Lindblad resonance with an m = 4 pattern from the 
effects of an ultra-harmonic resonance with an m = 2 pattern. 
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Table 1. Parameters of the phase mixed DF, / m i x . 

Disc _R(j o>o u z Lo q 

(kpc) (kms -1 ) (kms -1 ) (kpckms -1 ) 



Thin 3.0 27 20 10 0.45 

Thick 3.5 48 44 10 0.45 
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Figure 1. Figures showing the lines in action space associated 
with an 2:1 ILR (dotted), 4:1 ILR (solid), 2:1 OLR (short- dashed), 
and a 4:1 OLR (long-dashed) in the Galactic potential used in this 
study (upper) and a logarithmic potential (lower), for perturbcr 
pattern speeds Q p chosen such that the resonance includes the 
point J r = 0, J,/, = 2.09 (note that in all cases except the 2:1 ILR 
this is the resonant line used in the Hyades model) . In each figure 
the distribution in action of stars observed by the GCS is shown 
as a contour plot, both to give a sense of scale and to show the 
Hyades in action space (the broad overdensity near the resonance 
lines at higher J r ). Changed Q p for any of the resonances has the 
effect of moving the resonance line in while having a negligible 
effect on the gradient of the line. This figure illustrates the fact 
that the resonant lines in action space are all very similar to one 
another, and can change significantly with a change of assumed 
potential. 



(and therefore a flat rotation curve) and assuming a typi- 
cal Galactic potential (taken from McMillan 2011a) which 
has a circular speed which varies by ~ 4kms~ 1 in the 2 kpc 
either side of the Solar radius. It may at some point be possi- 
ble to use these differences to apply observations of resonant 
structures to accurately determine the shape of the Galactic 
potential (Binney 2005), but that is well beyond the scope 
of this study. 

It is also clear that the stars associated with the Hyades 
are not uniformly distributed in angle - most clearly because 
they are all moving radially outwards in the disk whereas, 
if they were uniform in angle, there would be as many stars 
moving out as in. This can also be explained in terms of 
resonant trapping (e.g. Binney & Tremaine 2008, section 
3.7.2). If we define 6$ — 9$ — Q p t (effectively the <j> angle 
coordinate in the frame rotating with pattern speed fl p ), 
then it is clear from eq. 1 that for orbits near resonance, 
the value of W r + mB<f, will evolve slowly - this is known 
as the "slow angle" - and the dynamical effect of the non- 
axisymmetric perturbation is dominated by forces that affect 
only the slow angle (as other forces rapidly average to zero) . 
Resonantly trapped orbits are ones in which the slow angle 



librates about a fixed point, with an amplitude that can be 
of order unity (Binney & Tremaine 2008). 

Because of this effect, Sellwood (2010) pointed out that 
one should also expect to observe an overdensity in angle 
for any set of stars trapped at resonance about the line in 
angle space 

W r + mOtf, = const, (2) 

where I and m take the same values as in eq. 1. Sellwood 
argued that the Hyades, as seen by the GCS, was associated 
with an overdensity about one of these lines with I = — 1 
and m — 2, 3 or 4, which would be a clear indication that it 
is a inner Lindblad (or ultra-harmonic) resonance. 

However, Mil showed that the selection effects associ- 
ated with sampling a finite volume (such as the Solar neigh- 
bourhood as observed by the GCS) play a complicated role 
in shaping the observed distribution in angle space. This 
means it is too difficult to disentangle the true signal asso- 
ciated with some distribution in described by eq. 2 (from 
which one might determine the values I & m) from the signal 
associated with the distribution in J combined with the se- 
lection effects. This conclusion was supported by Hahn, Sell- 
wood, & Pryor (2011) who investigated the Hyades in the 
Solar neighbourhood (within 200 pc) as observed by RAVE 
and the Sloan Digital Sky Survey (Abazajian et al. 2009). 

None the less, it is clear that the different distributions 
in angle associated with eq. 2 are the key to discriminat- 
ing between the different resonances that may be associ- 
ated with the Hyades, as the resonances are associated with 
nearly identical areas in action space. It is, however, neces- 
sary to look beyond the immediate Solar neighbourhood. 



3 NUMERICAL DETAILS 

In all cases we use the "convenient" model Galactic poten- 
tial given by McMillan (2011a). This model consists of a 
bulge component, thin and thick exponential discs, and a 
Navarro, Frenk, & White (1996) halo. This sets the solar ra- 
dius Ro — 8.5 kpc and the circular velocity at the Sun (the 
local standard of rest) vo = 244.5 km s -1 . When we need 
to transform observational data into Galactocentric velocity 
measurements we assume that the velocity of the Sun with 
respect to the local standard of rest is the best-fitting value 
found by Schonrich, Binney, & Dehnen (2010) 

do = (U@,V Q ,W e ) = (11.1, 12.24, 7.25) kms"\ (3) 

with U being the velocity towards the Galactic Centre, V 
being the velocity in the direction of Galactic rotation, and 
W the velocity towards the north Galactic pole. 

We use the torus-fitting method to find the values for 
x, v for any values of 0, J. All actions are quoted in units 
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Table 2. Parameters of the resonant component of the DF, / res - These are the parameters used in eqs. 8, 9 and 10. Q p and -Rcorotation 
are the pattern speed and corotation radius of the perturbation modelled in each case. Actions are in units of kpc 2 Myr -1 . 



Resonance 


^(Myr- 1 ) 


-^corotation ( kpc) 


J<t> 


,res 




B 


jR,cut 


^R,cut 


Co 




OLR 2:1 


0.058 


4.5 


2.09 - 


3Jr/2 


0.05 


0.3 


0.004 


0.001 


-1.7 


0.3 


OLR 4:1 


0.051 


6.0 


2.09 - 


SJr/5 


0.05 


0.3 


0.004 


0.001 


-2.0 


0.3 


ILR 2:1 


0.028 


28 


2.11 


-3J r 


0.05 


0.3 


0.004 


0.001 


1.3 


0.3 


ILR 4:1 


0.036 


13 


2.09 - 


9J R /5 


0.05 


0.3 


0.004 


0.001 


1.0 


0.3 



of kpc 2 Myr -1 . 4 The zero-points for the angles can be de- 
fined arbitrarily (provide the same convention is applied for 
all orbits), and for clarity we follow the conventions used by 
Mil, so each component of 9i lies in the range [— 7T, 7t], we 
define the zero point of 9$ such that at apocentre 9^ — (j>. 
and we take 9 r — at apocentre, and therefore 9 r = ±7r at 
pericentre. The Galactocentric coordinates are aligned such 
that the Sun is at a position in real space with Galactocen- 
tric coordinate = 0. Note that for small J r , J z , where it is 
appropriate to use the epicycle approximation, the value of 
9$ corresponds to the position (in <fi) of the guiding centre. 

The dynamical models we use have two elements - a 
phase mixed component with a DF /mix (J), and a compo- 
nent that is associated with the resonance which produces 
the Hyades with df f Tes (0,J). Both are based around the 
"pseudo-isothermal" dfs which have been used in a number 
of recent papers to describe disc dynamics in the Milky Way 
(Binney 2010; Binney & McMillan 2011; Binney 2012, Mil). 
These are of the form 



/(Jr, J<t>,Jz) = fa r (Jr, Jtfr) X 



27ral 



where 



fa r ( Jr i J <$> ) — 



[1 + tanh(J^/L )]e" 



(4) 



(5) 



R c 



Here k(J$) and v(J<j,) are the radial and vertical epicycle 
frequencies respectively and E(J^) = Y. e~ < - Rc ~ R °^ Rd is the 
(approximate) radial surface-density profile, where Rc(J<j>) 
is the radius of the circular orbit with angular momentum 
J^>. The factor 1 + tanh(J^/Lo) in equation (5) is there to 
effectively eliminate stars on counter-rotating orbits and the 
value of Lo is unimportant for these examples provided it is 
small compared to the angular momentum of circular orbits 
at the radii of interest. a r and a z are both functions of 
which control the radial and vertical velocity dispersions, 
with 



&z(J<j>) 



cr z o e 



q(R -R <: )/R d 



(6) 



where cr r o,cr z o and q arc constants. 

The phase mixed component with DF / m ix(J) is iden- 
tical to that used by Mil. It is the sum of two "pseudo- 
isothermal" dfs, corresponding to the thin and thick discs, 
normalised such that at the Sun the surface density of thick- 
disc stars is 23 per cent of the total stellar surface density. 
Table 1 lists the parameters of each component of the DF. 



4 N.B. 1 kpc 2 Myr 1 978 kms 1 kpc. The angular momentum 
of a circular orbit at the Solar radius is ~ 2.13 kpc Myr - 1 = 



The resonant component of the DF, f rcs (G,J), is of a 
similar form to that used by Mil, and can be written as 



/ rcs (6>, J) = A rcs / mix (J) x a(J) x /3(e) 



(7) 



A rcB is a normalisation constant, chosen such that in the 
Solar neighbourhood (taken to be a sphere of radius 200 pc 
about the Sun), 8 per cent of the stars are associated with 
the resonance. The term a(J) describes how the stars are 
trapped around the resonance line in action space. Sellwood 
(2012) shows that for particles at an ILR in an idealised N- 
body disc this trapping leads to an excess of particles near 
the resonant line at high Jr, and a dearth of particles at low 
Jr. We assume that this is what we would expect for stars 
affected by any of the resonances considered. This leads us 
to use 



t I {'Jtfi J<p^cs(Jr)) 

q(J) = exp ( -±-2— — \ K —LL_ | 

J, res 



(8) 



-B + -(l + B)(l + tanh(( Jjj - J R , C ut) / Aj, cut )) 



which is a Gaussian in of width Aj re3 centred on 
J<t>,res(Jr), multiplied by a function which runs from 1 at 
Jr > JR,cut to — B at Jr, < jR, cu t with the transition be- 
ing over a range controlled by the parameter Aj ;Cut . The 
line = J<f,,res(Jr) is, in each case, a good approximation 
to the line described by the condition on orbital frequency 
(equation 1) for a given perturber pattern speed O p . We are 
careful to ensure that in all models / m ix(J) + fies(J) > 
for all J. 

The term f}(&) in equation 7 is 



13(0) = exp 



(9r 



A 2 

^-V res 



Jt/>,res(Jr) is chosen such that 
m£lj > (J ri J^.res) = const for J z = 
and is chosen such that 

I9r,res + m9 t f, = Cg, 



(9) 

0, and 9 r ,res(9<l > ) 

(10) 



2080 kms' 1 kpc. 



where the constant Ce = I9 r ,res(94, = 0). The values Aj, res 
and Ag^es give the width of the resonance peak around the 
exact resonance lines in J^ and 9 r , respectively. 

The parameter values for the resonances considered in 
this study are chosen to provide a reasonable match (by eye) 
to the Hyades as seen in the GCS, and are shown in Table 2. 
This matching is a relatively simple process as the approx- 
imate positions of the resonance lines in action space are 
easy to determine from the density distribution of the GCS 
stars in action (Figure 1), and the approximate positions of 
the resonance lines in angle (and thus the value of Ce) from 
the density distribution in 9 r — 9$ plane (or marginalised 
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Figure 2. Contour plots of the density in the (— vr) — vj, plane 
of the stars observed by the GCS (black) and stars in the 200 pc 
around the sun in each of the models described in Table 2 (red) . 
— vr is used so that the plot reflects the familiar U — V plane 
in the Solar neighbourhood. The Hyades are found around v r = 
30kms -1 , ?)0 = 240kms _1 , and each of the models matches the 
observed overdensity well. 



Figure 3. Distribution of stars associated with the Hyades in 
the plane of the Galaxy for each of the four models (in each case 
these are 5000 stars selected at random from the positive part of 
/res)- The arrow indicates the direction of Galactic rotation. The 
dashed line is drawn at the Solar radius. The nine bins on the left 
hand side of each diagram are the bins shown in Figure 4, so the 
Sun is at the middle of the central bin, and each bin spans 600 pc 
in radius and 6.1° in azimuth. 



versions thereof, e.g. Fig. 8 or 5 of Mil). Contour plots of 
the density of stars in the (—vr) — Vcf> plane in the Solar 
neighbourhood in the four models, overlaid on the same fig- 
ure plotted using the GCS data are shown in Figure 2. This 
emphasises the point made by Mil that the information 
provided by the distribution of stars in velocity space in the 
Solar neighbourhood is insufficient to determine the type 
of resonance which gives rise to the Hyades moving group. 
It is also worth noting that very similar local velocity dia- 
grams can be produced with models that have significantly 
different values for various parameters (e.g. either Aj jres or 
Ae.res can vary by a factor of ~ 2. 



4 BEYOND THE SOLAR NEIGHBOURHOOD 

It is clear that to determine the resonance responsible for 
the Hyades moving group, we need to look beyond the Solar 
neighbourhood. In this section we determine what each of 
the Hyades models predicts for the distribution of stars in 
velocity space at different points in the Galaxy. 

In Figure 3 we show the positions of stars associated 
with the Hyades in the Galactic plane for all models. Much 
as one would expect from the condition described by equa- 
tion 2, the 4:1 resonances have groups of stars gathered at 
apocentre or pericentre at four different azimuths (8 r chang- 
ing by 8n as 6$ changes by 2ir), producing a square or clover- 
leaf pattern in position in the Galactic plane. Similarly the 



5 It is more usual to show this plot in U and V but, since we also 
consider fields far from the Solar neighbourhood, we use — vr and 
which, at the Sun's position, are equivalent to U and V except 
offset by the Sun's velocity 



a 
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Figure 4. Spatial bins in R & used in the plots shown in 
Figures 5 and 6 (solid lines). Dotted lines run in the R and (f> 
directions. The Sun is marked (as 0) at R = 8.5, <f> = 0. 



2:1 resonances have groups at apo- or peri-centre at two dif- 
ferent azimuths (8 r changing by 4-7T as 6$ changes by 2-7r), 
producing an oval pattern in position. 

In Figure 3 we also show the position of the extended 
Solar neighbourhood which we will be investigating in ve- 
locity space to find the signature of the Hyades. The key 
thing to note is that the "stream" in position associated 
with the resonances runs in different directions near to the 
Sun, depending on whether they are inner or outer Lind- 
blad resonances. ILRs produce streams that tend to be at 
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R < Ro in the immediate anti-rotation direction, and to- 
wards R > Ro in the direction of rotation, whereas OLRs 
produce streams that are the other way around. Note that 
each case the Hyades stars in the Solar neighbourhood are 
all rotating in the same direction as the Galaxy and moving 
radially outwards, so they are not "following" the stream 
in the Galactocentric frame (though they are in the frame 
rotating with the perturbation). 

These streams are of non-negligible width, controlled by 
the parameter Ag tleB in equation 9. It should be possible to 
determine the value of Ao, r es (or something similar, for a 
different model), as well the type of resonance, from obser- 
vations by determining the range in R (at a given <j>) over 
which the Hyades can be found (in velocity space). Sellwood 
(2010) attempted to estimate the time since the resonance 
occurred by comparing the width of the resonance in fre- 
quency to the width of the resonance in angle. That result 
is flawed because of the unrecognised impact of selection ef- 
fects on the observed spread in angle, but determining the 
true value of A.e,res would enable us to make an improved 
estimate. It should be noted that this can only be an upper 
estimate, without knowing the width in intrinsic in the 
libration of the resonantly trapped stars. 

In Figure 5 we show the distribution in velocity space of 
stars associated with the Hyades in each model, separated 
into bins in the Galactic plane each of which covers 600 pc 
in Galactocentric radius R, 6.1° in Galactocentric angle <f> 
(corresponding to ~ 900 pc at the Solar radius), and only 
contains stars which are within 300 pc of the Galactic plane. 
These bins are illustrated in Figure 4. 

In Figure 5 the Hyades can clearly be seen as a strong 
feature in all three azimuthal bins at the Solar radius for all 
models. In each case the Hyades appear at a similar point 
in the v r — v^ plane as they do in the Solar neighbourhood. 
This is in keeping with the overdensity observed by A12. In 
bins at lower R, the Hyades are found at higher vj,, and in 
bins at higher 7?, they are found at lower v$. This is be- 
cause in all cases the Hyades are associated with a narrow 
range in angular momentum, which naturally produces a 
correlation between R and v$. The peak of the distribution 
associated with the Hyades moves from v$ ~ 250 km s -1 in 
the inner bin to ~ 230 km s -1 (varying slightly between 
the models). This is a somewhat smaller change than one 
would naively expect from angular momentum conservation 
given that the bin centres are 600 pc apart, but this is be- 
cause stars in the outer bins tend to be close to the central 
bin rather than spread evenly across the bin. 

There are two notable differences between the ILR mod- 
els and the OLR models. Probably the most useful differ- 
ence, in terms of telling one model from another, is the vari- 
ation in the number of stars associated with the Hyades in 
each spatial bin. The two ILR models tend to have many 
Hyades stars in the outer Galaxy, and few in the inner, 
in the direction of Galactic rotation and many in the in- 
ner Galaxy, and few in the outer, in the anti-rotation di- 
rection. The opposite is true in the case of the OLR mod- 
els. This effect is most pronounced in the 4:1 resonances, 
with the 2:1 resonances placing few Hyades stars far from 
the Solar radius in this range of <j>. This is because of the 
differences in the paths of the "streams" illustrated in Fig- 
ure 3. The second difference is in the v T value of the peak 
in the velocity distribution in the inner or outer bins. The 



ILR models have a peak at v r ~ 50kms _1 in the bins 
7.6 kpc < R < 8.2 kpc and v r ~ 22kms~ 1 in the bins 
8.8 kpc < R < 9.4 kpc, whereas in the OLR models this 
is reversed, with v r ~ 22kms~ 1 for 7.6 kpc < R < 8.2 kpc 
and v r ~ 50kms _1 for 8.8 kpc < R < 9.4 kpc. 

In Figure 6 we show the same figures for the model 
including the smooth background df. This shows that in 
the inner and outer bins with relatively strong signatures of 
the Hyades, we can expect to observe the Hyades above the 
background, but the weak signatures are lost in the back- 
ground. It should be noted that this is, in some ways, a 
best case scenario, as we have assumed negligible uncer- 
tainty on velocity, a large number of observed stars (so shot 
noise doesn't hide the Hyades) and ignore other substruc- 
ture in the velocity diagram that might cause confusion (for 
example, from the dynamical structures associated with the 
Pleiades and Sirius moving groups). 

We can look more closely at the distribution of stars 
in spatial bins similar to those used by A12. If we take a 
Galactocentric Cartesian coordinate system (X,Y,Z), with 
the Sun at X — Ro — 8.5 kpc, Y = 0, and with Y increasing 
in the direction of Galactic rotation at the Sun, then we can 
take stars with —1.3 kpc < Y < —0.7 kpc, and in windows in 
the range 7.6 kpc < X < 9.4 kpc to correspond to the areas 
probed by A12's work (though they chose to take Ro = 
7.8 kpc). To probe the area in detail, while still using large 
spatial bins to maximise the numbers of stars per bin in 
a given survey, we consider overlapping bins. Each bin is 
0.6 kpc wide in X, and we consider bins with limits that 
are shifted by 0.2 kpc in X from each other. These bins are 
illustrated in Figure 7. 

In Figure 8 we show the distribution of stars in the 
(—vr) — space in each of these spatial bins. Again it is 
clear that the major difference between the various models 
is whereabouts in the Galaxy the Hyades remains a strong 
feature. Since these bins are in the anti-rotation direction 
from the Sun, the OLR models show a strong Hyades feature 
in the outer bins, and the ILR models show it in the inner 
bins. 

Finally we look at the consequences of observational 
uncertainties on the structures observed in the (—vr) — v^, 
plane for the various models. If we assume an uncertainty of 
2kms _1 in line-of-sight velocity, lmasyr -1 in proper mo- 
tion, and 20 per cent in distance (approximately the uncer- 
tainty on RAVE observations), and apply this to our models, 
we find the "observed" distribution in the (— vr) — plane 
shown in Figure 9. The most obvious consequence is that 
the entire velocity distribution gets stretched in the direc- 
tion associated with the uncertainty in transverse velocity. 
This is because the dominant velocity uncertainty is that 
associated with the distance uncertainty. This uncertainty 
also weakens the signature of the Hyades in each case. The 
differences between the two models is discernible, but it is 
seriously faded. It is worth noting that the signatures of 
the ILR models are visible in the innermost spatial bin, but 
those of the OLR models are not in the outermost. 

4.1 Comparison to observations 

Isolating these signatures in real data is likely to be even 
harder than Figures 8& 9 would suggest, because there will 
be additional noise due to the relatively small number of 
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Figure 5. The four different panels correspond to the four Hyades models (as labelled). Each panel is divided into nine plots, correspond- 
ing to the nine spatial bins shown in Figure 4, and shown in the same orientation (i.e. plots in the left of each panel all correspond to bins 
between 8.8 and 9.4 kpc from the Galactic centre and those at the bottom of each panel all correspond to the range —9.15 < <t> < —3.05). 
Each plot is a contour plot of the distribution of stars associated with the Hyades (i.e. / res and not f m i x ) in the (— vr) — plane. In 



each individual plot the range in — vr is —60 to 60kms — 1 and that in is 180 to 300 kms -1 (i.e. the same range as Figure 2), with 
dotted lines marking vr = Okms -1 and = 244.5kms — 1 = vo, the circular speed at the Sun. Contours are placed in each plot where 
the density is 10, 30, 50, 70 or 90 per cent of the highest density seen in the central plot of the panel (the bin which contains the Solar 
neighbourhood). Red crosses mark the highest density point in each plot. 



stars that will have been observed in any given volume, and 
because the other stars will not be smoothly distributed in 
velocity space, but will be in various substructures, much 
as they are in the Solar neighbourhood. However, sophisti- 
cated analysis techniques such as wavelet transforms (A12 
and references therein) may provide a means of finding even 
this weak signal in noisy, background dominated, data. 

Currently available data, as analysed by A12, does show 
a feature consistent with any of these models at radii near 
Ro- It is shown in their Fig. 3, second panel (labelled Sr q ), 
with the feature labelled "3" as it is the third greatest over- 
density in the plot. 

However, the panels showing the velocity distribution 
in bins at smaller or larger Galactocentric radii (labelled 
Sin and Sout respectively) do not show features that can be 
convincing associated with any of these models. The feature 
labeled "2" in the 5i n figure may appear to coincide with the 



overdensity expected from the ILR models shown here, but 
is at significantly higher v,p (~ 15 — 25 kms -1 faster than 
circular, as opposed to ~ 5 kms -1 for the ILR models), and 
lacks the extension towards the left of the plot (i.e. towards 
higher positive vr) that is a distinctive feature of all of these 
models (and the Hyades in the Solar neighbourhood). While 
there appears to be some small overdensity at the position 
that the OLR models would produce one in the S ou t figure, 
it is so minor that it is unlabelled, and is of completely the 
wrong shape (the feature labeled "4" in this figure is at much 
lower «0 than one would expect for the Hyades and, as noted 
by A12, is almost certainly the Hercules stream). 

A comparison of these models to figures showing the 
distribution of the overdensities in velocity space in more 
closely spaced, overlapping bins taken from the RAVE data 
(Antoja, priv. comm.) does not provide any significant fur- 
ther evidence for either the ILR or OLR models. Given the 
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Figure 6. The same plots as Figure 5 except that the full DF is used 
out the Hyadcs from a smooth background. 

substantial impact of observational errors illustrated by Fig- 
ure 9, and the further problems associated with low number 
statistics and a non-smooth background, this cannot reason- 
ably be taken as evidence that these models are all wrong, 
but rather as an illustration that more accurate data will 
be required to determine which, if any, of these models is 
correct. 



5 CONCLUSIONS 

In this paper we have taken a detailed look at four different 
models for the Hyades, each of which accurately reproduce 
the signature of the Hyades moving group in the Solar neigh- 
bourhood, but which differ significantly from one another 
beyond it. The models represent stars trapped at four dif- 
ferent resonances, and are produced using the torus fitting 
method to realize models described in angle-action variables. 

We have shown that for each model we expect the 
Hyades to produce a significant overdensity in velocity space 
around the Solar radius, even ~ 1.5 kpc from the Sun in 
Galactic azimuth. However this signature moves and fades 
with Galactocentric radius in a predictable way that pro- 
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(i.e. both fres and f m i x ), giving a sense of the difficulty in picking 

vides a clear way of determining which kind of resonance is 
responsible. Most notably, if the Hyades are due to an ILR, 
they tend to be found at larger Galactocentric radii (R) in 
the direction of Galactic rotation, and smaller R in the anti- 
rotation direction, whereas the opposite is true if they are 
due to an OLR. 

It is worth noting that, as shown by Sellwood (2012), 
there are good theoretical reasons to expect trapping at an 
ILR, and it can be seen in simulations. It is not obvious 
that one should expect to see similar at an OLR, and nei- 
ther trapping at the OLR or the ultra-harmonic resonance 
is observed in Sellwood's simulations. An additional concern 
is that, as noted by Sellwood (2010), a significant spiral pat- 
tern with co-rotation around 28 kpc, as would be required 
to produce the Hyades through a 2:1 ILR, seems unlikely, 
which means that it would have to be due to some other per- 
turbation, such as a rotating m = 2 distortion of the halo. 
If the Hyades are the result of a an inner ultra-harmonic 
resonance with an m = 2 perturbation (equivalent, in the 
models used in this paper, to a 4:1 ILR), this would be in 
keeping with the results of Siebert et al. (2012), who used 
model with a two arm spiral to explain a gradient in mean 
Galactocentric radial velocity found in RAVE data, and the 
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Figure 7. Spatial bins used in the plots show in Figures 8 and 9 (solid lines). The dotted lines run in the R and <j> directions at the 
same places as in Figure 4 to indicate scale. 




model proposed to explain some elements of the local ve- 
locity distribution by Quillen & Minchev (2005). It should 
be noted that the absence of resonant trapping at the ultra- 
harmonic resonance in the simulations of Sellwood (2012) 
would seem to argue that if this is the result of resonance 
with a spiral with corotation at ~ 13 kpc, the spiral would 
have to be 4-armed. 

Presently available data (A12) suffers from low number 
statistics and uncertain distances, which leads to a signifi- 
cant uncertainty in the transverse velocity of the stars. This 
means we are not yet able to provide any significant extra 
insight as to which resonance is associated with the Hyades 
using observations from beyond the Solar neighbourhood. 

With improved distance measurements and more ob- 
served stars, this situation will improve, and these easily 
tuneable models provided by the angle-action description of 
the phenomena will provide great insight into the nature of 
the resonance and, therefore, the perturbation to the Galaxy 
that is responsible. 
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